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Abstract 

We rederive the finite size scaling formula for the apparent critical temperature by 
using Mean Field Theory for the Ising Model above the upper critical dimension. We 
have also performed numerical simulations in five dimensions and our numerical data 
are in a good agreement with the Mean Field theoretical predictions, in particular, 
with the finite size exponent of the connected susceptibility and with the value of the 
Binder cumulant. 



1 Introduction 



In this letter we will study the scaling above the upper critical dimension for Ising Models (where Mean 
Field holds). This issue was been controversial because the numerical simulations do not support the Mean 
Field predictions 0|. In particular neither the scaling exponent of the apparent critical temperature, nor 
the scaling exponents of the connected susceptibility, nor the value of the Binder cumulant were those of 
the Mean Field [|[ ^ . The first discrepancy has been solved using the renormalization group D and we will 
show that it can also be explained using Mean Field Theory. We have also made numerical simulations 
in five dimensions and we see that our numerical values are in a good agreement with the Mean Field 
predictions and hence there is no discrepancy. 

The plan of this letter is the following. In the next section we will made some remarks on the different 
limits of the Binder cumulant. We will then derive the scaling formula for the apparent critical temperature 
using Mean Field Theory. Next we show the numerical results for the five dimensional Ising model and 
finally we will report the conclusions. 



2 Limits of the Binder Cumulant 

The problem of the triviality of the field theory can be understood studying the value of the renormalized 
coupling constant. We center the following discussion using one component spin models. 

If we renormalize the theory at zero momenta it can be shown that the renormalized constant, gn, 
obeys: 



'«ns:J^' (I) 

where ^oo is the correlation length defined through the decay of the correlation function, d is the dimension, 
and B, the Binder cumulant, defined as: 

where M is the total magnetization. With this definition if the distribution of M is Gaussian B vanishes. 

From this point of view, triviality is equivalent to have a renormalized coupling constant, gji, which 
vanishes in the thermodynamical limit at the critical point (where we define the field theory). But gn 
equal to zero implies that the Binder cumulant is also zero. This corresponds to requiring that the limiting 
(infinite volume) distribution of the total magnetization is Gaussian, as happens, for instance, in the high 
temperature phase. In a non trivial theory, the limiting probability distribution of the magnetization is 
not Gaussian, and the degree of the non-gaussianity is measured by the value of the renormalized constant. 
In the low temperature phase the limiting distribution is the sum of two Dirac deltas (one with positive 
magnetization and the another with negative) and then the Binder cumulant is one ^ . 

In this context is clear that we can use the value of g^ to classify the universality class and since L/^^l 
is a constant at (3c we could also use the Binder cumulant for the same purpose. 

The problem is that we measure the Binder cumulant on a lattice L at temperature 1/(3, and so we 
obtain a function B{L,(3). The finite size scaling behavior of the Binder cumulant is 

B{L,(3) = f{L^/%P^Pc)). (3) 

The scaling function has the following limits: 

• /? is fixed to the critical value (3c and then we take the infinite volume limit: 

Bi EE lim lim B{L,(3). (4) 
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• Wc take the value of the Binder cumulant at the apparent critical temperature, denoted as (3c{L) (de- 
fined, for example, as the maximum of the connected susceptibility, obviously /3c = liniL^oo /3c(i))j 
and we consider the limit: 

B2= lim B{L,I3,{L)). (5) 

L — >oo 



First take L — > oo and then /? — > /3c 



Bs = lim lim B{L,(3). (6) 

(i—^^c L—^oo 



In general, these limits are all different: 

Bi^B2^ B3. (7) 

Assuming that we have the following scaling for the apparent critical point, Pc{L) 

l3c{L)= + (8) 

where a is dependent on the precise definition of (3c{L) (for example the specific heat or the susceptibility). 
Then we have: 

B2 = f{a). 

Obviously B2 depends on the observable used to define the apparent critical temperature. This limit can 
be related to the usual calculations of the Binder cumulant at L/^l — const: at the maximum of the 
connected susceptibility in a finite size we have that the ratio between and L is constant, independent 
of the size. 

We remark that if we rewrite equation (^ as 

B(L,/3)=/(L/Cl(/3)), (9) 

and calculate the Binder cumulant for a set of /3's and L's such that L/^l(/3) — const, the value of this 
Binder cumulant will be independent of the /3's and L's in the scaling region (defined as the region where 
the formula (^, or equivalently (^, holds). 
For completeness: 

Bi =/(0) , i?3-/(oo). (10) 

Bi and B2 are both size independent and observable independent: are the zero value and the infinite value, 
respectively, of a universal function. 

In next section we will show how B2 can be calculated using the connected susceptibility to define 
/3c(L). 

3 Above the critical dimension. 

We will follow the notation of reference in our presentation. Above the critical dimension we can study 
the system considering an effective action with only one degree of freedom 

5eff = L^' + , (11) 

where t is the reduced temperature and we need to calculate: 

iyd</>02fexp(-^eff), (12) 



m2p ^ 



where 



Z = J d<^exp(-S'eff). (13) 



3 



In this notation the Binder cumulant is simply 

(3-TO4/to^)/2. 



(14) 



We perform the following change of variables: 

$ = (uL'*)-i/V, (15) 

and the action becomes 

S,f,^^xit,L)<^' + ^<^\ (16) 

where x{t,L) = ti''/2y-i/2^ 

We remark that at this moment we are in a finite size and at a generic temperature, different from the 
critical one. 

If we calculate the value of x{t, L) at the maximum of the connected susceptibility {< (fp' > ~ < |(/)| >^) 
with the action of the equation ( p^ ) (in other words, we calculate the apparent critical reduced temperature 
t{L)) we obtain, evaluating the integrals numerically: 

x{t, i)max = -0.44015 = t{L)L''^^u-^/^. (17) 

and hence the scaling of the apparent reduced temperature is 

t{L) - L-'^/^ (18) 

as was previously found in reference Q| and numerically in ||l|. The subdominant term L^^'' follows from 
the subsequent analysis of Brezin and Zinn- Justin Q. The analysis of Brezin and Zinn- Justin was done 
with t = and they miss the dominant term L^'^/'^ because it is proportional to t. 

We can evaluate the Binder cumulant, eq. (|lj), at this apparent critical temperature with the action 
( prj ) obtaining 

^2 = S(L,/3c(i)) = 0.742137. (19) 

This value is independent of the dimension, size and of the coupling u but depends of the observable used 
(in this case the connected susceptibility). The generalization to 0{N) models is straightforward. 



4 Numerical Results in = 5. 

We have performed numerical simulations in order to check the discrepancy between previous quoted values 
and the theoretical prediction. 

To do this we have simulated the five dimensional Ising Model using the Wolff single cluster algorithm. 
We have simulated L = 4, 6, 8, 10, 12 and 16 for a total number of cluster sweeps of the order 6.5 million. 
In the analysis of the data we have used the spectral density method |0] and the jack-knife method to 
evaluate the statistical error. 

The shift of the critical temperature, calculated as the maximum of the connected susceptibility, follows 
(using all the sizes in the fit) 

l3c{L) = 0.11387(4) + 0.0798(20)L"2.32(i6)^ (20) 

with a x^/d-o.f « 0.2 and we write the errors inside of the parenthesis. Hence we obtain a shift compatible 
with the theoretical prediction 5/2. This value was also reported in ||l|. We have repeated the global fit 
with L >Q the numbers are very similar {(3c — 0.11386(6)) but the errors are large. 

Fixing the scaling exponent to the theoretical value of eq. (2.5) and using L > 8 , we obtain: 

l3c{L) = 0.11388(3) + 0.12(l)L-2-^ (21) 

with a good x^/d.o.f. « 0.15. In Fig. 1 we plot the numerical values of (3c{IA and plot, signed by a 
continuous line, the fit (|2|). We remark that the estimation of Rickwardt et al was f3c = 0.113929(45) 
that is a 1.5 standard deviations of our value. 
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Figure 1: Effective inverse temperature, (3c{L), against L for six different sizes. The continuous line is the 
global fit supposing that the exponent is 2.5 (as described in the text). 



L 




B{L,^,{L)) 


4 


0.485(5) 


0.739(3) 


6 


0.475(8) 


0.744(4) 


8 


0.450(2) 


0.739(7) 


10 


0.425(5) 


0.746(4) 


12 


0.39(2) 


0.73(1) 


16 


0.40(2) 


0.74(1) 



Table 1: B{L,(3c) and B{L,(3c{L)) for different sizes. The Mean Field theoretical value is 0.40578 for the 
first Binder cumulant and 0.742 for the second one. 
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Figure 2: Scaling of the connected susceptibility for six different sizes. The continuous line is the fit using 
only the points L > 8, as described in the text. 




Figure 3: Wc plot the curves of the Binder cumulant for L = 10 (•), L = 12 (*) and L = 16 (o). The 
horizontal line is the Mean Field prediction (B=0.406). The vertical line is our estimation of the inverse 
critical temperature. 
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Using /3c — 0.11388 we can monitor B{L, Pdoo)) that we plot in Table 1. For the large lattices the 
value is compatible with that reported by Brezin and Zinn- Justin in reference S = 0.40578. Also we 
have written in this Table the values of B(L, (3c{L) and the agreement is again very good with our previous 
calculated value (p^. 

Analyzing the maximum value of the connected susceptibility we found the following power law using 
all the lattices in the fit: 

Xm..-L^-'''^'\ (22) 
but the fit is very poor (x^/d.o.f. = 7). However fitting only the lattices with L > 8 we found 

Xr...-L'-''^'\ (23) 

with x^/d.o.f. w 1.95/2, and so the agreement is very good with the theoretical prediction 5/2. We plot 
the numerical data with the fit ( p3| ) in Fig. 2. 

We finally report in Fig. 3 the different curves of the Binder cumulant near the critical temperature. 
Using the critical temperature of reference we have found, using our L = 16 data, a value of the Binder 
cumulant B = 0.52(2) that is a two standard deviations of B = 0.479(25), the value reported by this group 

0- 

5 Conclusions. 

We have derived using Mean Field the finite size shift of the apparent critical temperature for the Ising 
model above the upper critical dimension. We also have shown that the numerical data for the five 
dimensional Ising model are compatible with the Mean Field predictions, in particular, the value of the 
Binder cumulant at criticality, the finite size exponent of the connected susceptibility and the shift of the 
apparent critical temperature agree very well with Mean Field Theory. 
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